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For the first time tlie experiments ALEPH, DELPHI, L3 and OPAL have presented preliminary 
results for fermion-pair production in e^e~ collisions on the full data set above the Z pole. A 
combined analysis of the Bhabha scattering measurements is performed to search for effects of 
contact interactions. In the case of two axial-vector (AA) currents the best fit to the data is 2.6 
standard deviations away from the Standard Model expectation, corresponding to an energy scale 
A — lO.Sli g TeV for contact interactions. For other models no statistically significant deviations 
■ are observed, and the data are used to set lower limits at 95 % confidence level on the contact 

^— «j ' interaction scales ranging from 8.2 to 21.3 TeV, depending on the helicity structure. 
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■ I. INTRODUCTION 



The Standard Model (SM) is very successful in confronting the data coming from the highest energy accelerators. 
Still, there are theoretical reasons to expect that it is an effective theory, valid in a limited energy range, and one 
1^ of the first questions in the quest for physics beyond the Standard Model is what is the relevant scale, where new 
\^ , phenomena will give experimental signatures, fn this paper we will follow a data-driven approach and analyze the 
y—^ full set of measurements on the reaction 
^ . 

O. e+e-^e+e-(7) (1.1) 

^ ■ 

' at centre-of-mass energies ^/s above the Z resonance from 183 up to 208 GeV. Rather than concentrating on a 
■ particular model, we will look for something unexpected. Bhabha scattering is chosen as a very sensitive probe, 
^ i' affected in many new physics scenarios. This analysis is a continuation of the work presented in [Q, based on the 
JL published differential cross sections for Bhabha scattering at energies 183 and 189 GeV. 

^ In the Standard Model the production of a fermion-pair // in e^e~ collisions is described by the exchange of 7 or 

(~| Z in the s-channel, and if the final state is identical to the initial one, also in the i-channel. The interest in studying 
fermion-pair final states above the Z pole at LEP2 is driven by the fact that many types of new physics scenarios 
. !^ ' can contribute to these processes. For this to happen, the couplings to the initial and final states should be different 
from zero. In the case of Bhabha scattering we just need a coupling to the electron and the positron. Even if the 
5-H Standard Model extension operates at an energy scale much higher than the accessible centre-of-mass energy for a 
. 5^ I direct observation, it can still give measurable effects by modifying the differential cross section through interference 
with the SM amplitudes. 



II. CONTACT INTERACTIONS 



Contact interactions offer a general framework for describing a new interaction with typical energy scale A 3> \fs. 
The presence of operators with canonical dimension > 4 in the Lagrangian gives rise to effects ~ 1/A^~^. Such 
interactions can occur for instance, if the SM particles are composite, or when new heavy particles are exchanged. 

For fermion-pair production, the lowest order flavor-diagonal and helicity-conserving operators have dimension 
six 1^. The differential cross section takes the form 

^ = SM{s,t)+e-Cint{s,t) + e^ -Cciis^t) (2.1) 
ail 

where the first term is the Standard Model contribution, the second comes from interference between the SM and the 

contact interaction, and the third is the pure contact interaction effect. The Mandelstam variables are denoted as s, 

2 

t and u. Usually the coupling is fixed, = 1, and the structure of the interaction is parametrized by coefficients for 
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16 1 collaborations have presented results for the total 



the helicity amplitudes: \riij\ < 1, where = L, R) labels the helicity of the incoming and outgoing fermions. We 
define 

where the sign of rj enables to study both the cases of positive and negative interference. 

III. EXPERIMENTAL DATA 

Measurements on fermion-pair production for the full data set of the LEP2 collider have become available recently. 
They include preliminary data from 192 to 208 GeV centre-of-mass energies and published results at 183 and 189 
GeV. In the following we will concentrate on the measurements of Bhabha scattering, where large data samples have 
been accumulated during the very successful LEP runs from 1997 to 2000. 

The ALEPH [|-|, DELPHI [|-|, L3 0-[l|] and OPAL [[3|- 
cross section and the forward-backward asymmetry ApB oi Bhabha scattering. In the case of DELPHI, L3 and OPAL 
the results are for all energy points and the scattering angle 6 is the angle between the incoming and the outgoing 
electrons in the laboratory frame. In the ALEPH case the forward-backward asymmetry is available only at 183 GeV 
and the scattering angle 9* is defined in the outgoing e''~e~ rest frame. The acceptance is given by the angular range 
-0.9 < cos 6** < 0.7 for ALEPH, by U° < 6 < 136° for DELPHI and L3, and by | cos6l| < 0.7 for OPAL. 

The experiments use different strategies to isolate the high energy sample, where the interactions take place at 
energies close to the full available centre-of-mass energy. This sample is the main search field for new physics 0. 
DELPHI, L3 and OPAL apply an acoUinearity cut of 20°, 25° and 10° respectively. ALEPH defines the effective 
energy, •\/s', as the invariant mass of the outgoing fermion pair. It is determined from the angles of the outgoing 
fermions. For details of the selection procedures, the statistical and systematic errors we refer the reader to the 
publications of the LEP experiments. 

IV. ANALYSIS METHOD 

The Standard Model predictions for Bhabha scattering are computed with the Monte Carlo generator BHWIDE [ p^ . 
We assign a theory uncertainty of 1.5 % on the absolute scale of the predictions |2^]. In all cases the individual 
experimental cuts of the selection procedures and the isolation of the high energy samples are taken into account. The 
results are cross-checked with the semi-analytic program TOPAZO ||2^ and the Monte Carlo generator LABSMC . 

The effects of new phenomena are computed as a function of the parameter e, defined in Section |l|. Initial- 
state radiation (ISR) changes the effective centre-of-mass energy in a large fraction of the observed events. We take 
these effects into account by computing the first order exponentiated cross section following |^3|. Other QED and 
electroweak corrections give smaller effects and are neglected. 

In total we have 57 data points: 32 from the cross sections (eight energy points and four experiments) and 25 
from the forward-backward asymmetries. A fitting procedure similar to the one in ||2^ , ^5|| is applied. A negative 
log-likelihood function is constructed by combining all data points at the eight centre-of-mass energies: 



, , / [predictionfSM, e) — measurement]^ \ 

-log^ = L( ^Ta^ ) (4-1) 

1-1 ^ r 

where prediction(SM, s) is the SM expectation for a given measurement (total cross section or forward-backward 
asymmetry) combined with the additional effect of contact interactions as a function of the scale A, measurement 
is the corresponding measured quantity and A — error[prediction(SM, e) — measurement]. The index r runs over all 
data points. The error on a deviation consists of three parts, which are combined in quadrature: a statistical error 
and a systematic error (as given by the experiments) and the theoretical error assigned above. The systematic errors 
account for small correlations between data points. The minimum of the negative log-likelihood function gives the 
central value of the fitted parameter e for each model, and the interval containing 68% of the total probability around 
the minimum is used to determine the values corresponding to one standard deviation in the positive and negative 
directions. 



^ For reviews see e.g. 
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V. RESULTS AND DISCUSSION 



The results of the fits for the different contact interaction models are summarized in Table |. In the same table 
the coefficients specifying the helicity structure of the investigated models are given. For all models the central value 
of the parameter e is no more than 1.05 standard deviations away from the Standard Model value e = 0, with one 
exception. The fitted value for the AA model is 

e = 0.0095 TeV-2 (5.1) 

or 2.6 standard deviations from the SM expectation. The quality of all fits is good, with values ~ 52/56 degrees 
of freedom. For the AA model the value is 46.8/56 degrees of freedom. The log-likelihood curves for the VV and 
AA models are shown in Fig. ^ 

In the cases where the data from the LEP collaborations show no statistically significant deviations from the SM 
predictions, we can derive one-sided lower limits on the scale A of contact interactions at 95% confidence level. This 
is done by integrating the log-likelihood functions in the physically allowed range of the parameters describing new 
physics phenomena, assuming a uniform prior distribution. The exact definition can be found in [ p4[ . The limits for 
positive or negative interference are summarized in Table |[ The results presented here improve on the limits obtained 



by individual LEP experiments P,n,pq,p7 14 1 or in a combined analysis W 



Now we turn to a discussion of the result for the AA model. Clearly this result is unexpected and requires careful 



analysis. The central value and the one standard deviation band of Eq. 5.1 correspond to a scale for axial- vector 
contact interactions 

A+ = 10.3 {'l:', ;;^;J TeV. (5.2) 

Let us use this central value in our investigation. We can reverse the view point and ask ourselves what are the 
expected deviations from the SM predictions for the total cross section and the forward-backward asymmetry of 
Bhabha scattering. The result as a function of the centre-of-mass energy is illustrated in Fig. |^. The AA model 
changes the total cross section by less than 0.5% in this energy range. The relative effect is practically constant for 
rising energy. For comparison, the combined error on the total cross section of the four experiments for the highest 
energy point (207 GeV) is 1.2%. An additional factor is the theoretical uncertainty of 1.5% on the absolute scale. 



discussed in Section IV. On the contrary, the forward-backward asymmetry is changed by (-0.0066) in absolute value 
(or -0.8% relative) at 183 GeV and by (-0.0083) or (-1.0%) at 207 GeV. The effect is rising with energy. It is clear that 
the AA model manifests itself mainly in the forward-backward asymmetry, which is not affected by the uncertainty on 
the absolute scale. Another favorable fact is that the statistical error on ApB is given by a{ApB) = ^(1 — Apg)/N, 
so for the same number of events N we have a reduction factor (~ 0.6) due to the large value of the asymmetry ~ 0.8. 
On the minus side this measurement requires recognition of the electron and positron charges, and hence stricter 
event selection asking for good tracks. Also corrections for charge confusion have to be applied. For comparison, the 
combined error on the forward-backward asymmetry for three experiments and the highest energy point (207 GeV) 
is 1.1%. As the ALEPH collaboration has not presented results on Aps above 183 GeV, the result discussed here is 
based on the asymmetry measurements of the other three collaborations. 

The measurements of Bhabha scattering from 192 to 208 GeV are preliminary. The systematic errors and the 
central values may still change. The charge confusion correction has the form 

A measured 

/^corrected FB (r o\ 

" l-2c ^^-^^ 

where c is the amount of events with wrong sign assignment. For example, if c = 0.05, it is enough to underestimate 
the charge confusion by 9% in order to underestimate the asymmetry by 1%. For c = 0.02, one has to underestimate 
the charge confusion already by 24% in order to underestimate the asymmetry by 1%. We take into account in our 
analysis this major source of systematic error, as given by the experiments. 

The forward-backward asymmetry can be affected by the interference between initial-state and final-state radiation, 
which is known to change the form of the differential cross section. This effect is estimated using the program TOPAZO. 
If we switch the interference on, Afb is increasing by 0.1% for an acoUinearity cut of 10° and by 0.07% for an 
acollinearity cut of 25°. So the effects are very small, and as the experiments use differential efficiency in the scattering 
angle, the impact on the measured asymmetry values is much smaller. We can conclude that the interference is under 
control. 

Another interesting point are the preliminary results on contact interactions presented by DELPHI , L3 j2^ and 
OPAL [n6[. The deviations of e from the SM value are positive in all cases and amount to 0.8, 0.7 and 1.4 standard 
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deviations respectively. All three collaborations observe a pull in the same direction, but their sensitivity is not enough 
to make any conclusion. In (l) a limit of A+ > 10.4 TeV at 95% confidence level is obtained, which covers about 40% 
of the favored region from this analysis at higher energies. 

In order to clarify the situation, the final results of the four LEP collaborations are needed. The best option is 
to combine the measurements of the differential cross sections in the four experiments, as illustrated in Fig. ^. The 
AA model gives a specific signature, enhancing the cross section in the central part and backwards, and only slightly 
reducing the forward peak. 

When interpreting the physical meaning of contact interaction scales, we should remember that a strong coupling 

2 2 

= 1 for the novel interactions is postulated by convention. If we assume a coupling of different strength, ^ — a' , 



the limits can be translated as A' — \fo! ■ A. For the extreme case of a coupling of electromagnetic strength the real 
scale can be around 1 TeV. 

The measurements of Bhabha scattering above the Z resonance reach already a higher level of precision than the 
best theoretical tools available. In order to fully exploit the physics potential of the large data samples collected 
during the LEP running from 1997 to 2000, improved theory predictions are very desirable. A combined effort of 
the four LEP collaborations is needed to answer the question: is Bhabha scattering opening the first window to new 
physics at the TeV scale? 
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FIG. 1. Log- likelihood curves for the VV and AA contact interaction models of the combined fits to the data on Bhabha 
scattering from the four LEP experiments. 
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FIG. 2. Effects of contact interactions on the total cross section and the forward-backward asymmetry for Bhabha scattering 
as a function of the centre-of-mass energy (upper plots) and difference to the Standard Model predictions (lower plots): SM 
- solid, VV model - dotted and AA model - dashed. The energy scale is A+ = 10.3 TeV in both cases. The angular range is 
44° < 6* < 136°. 
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FIG. 3. Effects of contact interactions on the differential cross section for Bhabha scattering for ceritre-of-rnass energy 206.8 
GeV (upper plot) and the ratio to the Standard Model predictions (lower plot): SM - solid, VV model - dotted and AA model - 
dashed. The energy scale is A+ = 10.3 TeV in both cases. The "data" points show the expected statistical error if we combine 
the measurements of the four LEP collaborations for the six energy points from 192 to 208 GeV. 

TABLE I. Results of contact interaction fits to Bhabha scattering. The numbers in brackets 
are the values of [r]LL,r]RRT'nLR,'nRL] defining to which helicity amplitudes the contact interaction 
contributes. The models cover the interference of contact terms with single as well as with a com- 
bination of helicity amplitudes. The one-sided 95% confidence level lower limits on the parameters 
A+ (A_) given in TeV correspond to the upper (lower) sign of the parameters r], respectively. 



Model 


Amplitudes 


e 


A_ 


A+ 




l'nLL,'nRR,'nLR,riRL] 


[TeV-2] 


[TeV] 


[TeV] 


LL 


[±1,0,0,0] 


0.0028 t°orott 


10.3 


8.3 


RR 


[0, ±1,0,0] 


0.0036 in^s 


10.2 


8.2 


LR 


[0,0, ±1,0] 


-0 0046 +0 0044 

-U.UU'lD -0.0046 


8.8 


12.7 


RL 


[0,0,0,±1] 


-0.0046 tr^tt 


8.8 


12.7 


VV 


[±1,±1,±1,±1] 


-0.0008 


18.0 


21.3 


AA 


[±l,±l,Tl,Tl] 


0.0095 inSfr 


16.5 


8.0 


LL-I-RR 


[±1, ±1,0,0] 


20 +0 0029 

U.UUZU _0.0035 


14.3 


11.9 


LR-I-RL 


[0,0,±1,±1] 


-0 0023 +0 0023 

-U.UUZO _0.0024 


12.4 


18.2 



O 

06 



1.1 
1.075 
1.05 
1.025 
1 

0.975 
0.95 
0.925 
0.9 



